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Abstract 

A new construction of complex Hadamard matrices of composite order d = pq, with 
primes p, q, is presented which is based on pairs of mutually unbiased bases containing 
only product states. We illustrate the method for many product dimensions d < 100 by 
analytically deriving complex Hadamard matrices, both with zero and non-zero defect. 
In particular, we obtain at least 12 new isolated Butson-type Hadamard matrices of 
orders ranging from 9 to 91. 

1 Introduction 

A square matrix H of order d is a complex Hadamard matrix if it is unitary, HH< = I, 
and if its elements have equal modulus. This definition generalises the concept of a real 
Hadamard matrix with matrix elements limited to the values ±l/y/d. The first known 
construction of such matrices is due to J.J. Sylvester [T] while they take their name from 
J. Hadamard who found that the determinant of a unitary matrix achieves its maximum 
if all its matrix elements have the same absolute value [2]. 

Since then, complex Hadamard matrices have made their appearance in various branches 
of both mathematics and physics. For example, they relate to the problems of finding 
bi-unitary sequences and cyclic n-roots [3], they can be useful in constructing certain 
*-subalgebras of finite von Neumann algebras and error correcting codes [5]. They 
also have applications in quantum information, representing an important ingredient in 
teleportation and dense coding schemes [6], and they are closely linked to mutually unbiased 
bases [7]- For a detailed overview of their applications, see [5JIH]- 

In view of their many uses, a complete classification of complex Hadamard matrices would 
be highly desirable but has not yet been achieved. All complex Hadamard matrices, up 
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to equivalence, are known for dimensions d < 5 [HIITO], but their classification remains 
incomplete for higher dimensions. Many continuous families of complex Hadamard matri- 
ces already exist for d = 6, including the three- and four-parameter families found in [IT] 
and [12], respectively; in addition, a single isolated complex Hadamard matrix has been 
found. 

General construction methods exist in composite dimensions |13[I14]. and continuous fam- 
ilies of complex Hadamard matrices have been obtained from so-called parameterisations 
of known Hadamard matrices )14H17| . There has also been some success in finding contin- 
uous families of complex Hadamard matrices for certain prime dimensions |18j . A survey 
of known complex Hadamard matrices is given in |19| for d < 16, with an updated online 
catalogue provided by |20| . 

In this paper we discuss a new technique to construct complex Hadamard matrices of order 
d = pq, p < q, with p and q prime, based on a link between Hadamard matrices and pairs 
of MU bases consisting entirely of product states. A (d X d) complex Hadamard matrix H 
is unitary, which is another way of saying that its columns represent orthonormal vectors 
of a basis in a Hilbert space of dimension d. By choosing its matrix elements to have 
modulus 1/yd, the columns of H will form a basis which is mutually unbiased (MU, for 
short) to the standard basis, since the squared modulus of the inner product of any vector 
from H with any vector from the standard basis equals 1/d. In other words, the matrices 
{/, H}, where / is the (d x d) identity matrix, represent a pair of MU bases. Since any 
pair of MU bases can be brought into standard form, which contains the standard basis, 
the classification problem of complex Hadamard matrices of order d is equivalent to the 
problem of finding all pairs of MU bases in the Hilbert space C d . 

The method we use to construct Hadamard matrices originates from previous studies of 
pairs of MU product bases in dimension six |21H23j . Generalising this method from d = 2x3 
to composite dimensions d = pq, we establish a general construction method resulting in 
previously undiscovered complex Hadamard matrices. 

The paper is organised as follows: In Section [2] we summarise those properties of complex 
Hadamard matrices which we will use in later sections. Section [3] includes our first main 
result, Theorem [3] which describes a general construction for complex Hadamard matrices 
of size d = pq. In Section |4] we apply Theorem |3] to find new complex Hadamard matrices 
for d < 15, and briefly touch upon higher dimensions as well as potential generalisations 
of the construction in Section |5] Section [6] contains a summary of our results. 



2 



2 Complex Hadamard matrices 



2.1 MU bases and Hadamard matrices 

Two orthonormal bases B = {|<fo)} and B 1 = in a Hilbert space C d are MU if and 

only if |(0j|-0,-)| 2 = 1/d for all i, j = 1 . . . d. In prime and prime power dimensions p n , 
n G N, complete sets of {p n + 1) MU bases exist )24[I25|. In composite dimensions such as 
d = 6, 10, 12..., the question of whether complete sets of (d + 1) MU bases exist remains 
open |26| . 

MU bases in a finite-dimensional Hilbert space C d are closely related to complex Hadamard 
matrices. Given a set of (r + 1) MU bases in standard form {I,B±, . . .B r }, with I the 
standard basis, the bases B±, . . . ,B r can be represented by (d x d) complex Hadamard ma- 
trices, Hi, . . . , H r . Here the vectors of each basis are given by the columns of the respective 
matrix. Since these matrices are MU to the identity matrix, their matrix elements have 
modulus l/\/d. 

Two complex Hadamard matrices H and K are equivalent to each other if one can write 
H = D1P1KP2D2, where D\,D2 are diagonal matrices with diagonal elements consisting 
of phase factors and P\,P2 are permutation matrices. As a consequence, every complex 
Hadamard matrix can be dephased such that the matrix elements of its first row and 
column equal \j\fd. 

A useful criterion to determine whether two complex Hadamard matrices are inequivalent 
is to calculate their Haagerup set, defined in |9]. 

Definition 1. The Haagerup set of a complex Hadamard matrix H of order n with matrix 
elements hij is given by 

A(H) = {hijhkih u hkj : i,j,k,l = 1, . . . ,n}, (1) 
where hij denotes the complex conjugation of hij . 

Two matrices with different Haagerup sets are inequivalent since the set A(H) is invariant 
under equivalences. Note, however, that two inequivalent matrices may have the same 
Haagerup set. 

Sets of MU bases are equivalent if they can be mapped to each other by any of the following 
unitary or anti-unitary equivalence transformations: a unitary map acting on all bases 
within a set; the multiplication of states within bases by arbitrary phase factors; the 
permutation of states within a basis, and the complex conjugation of all bases. In addition, 
one may exchange any two bases within a set. 

It is often the case that a single Hadamard matrix is contained in a subset of a larger 
family of complex Hadamard matrices of which there are two types: affine and non-affine 
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families. The following definition, given in 1 19] . uses the notation o as the entrywise product 
of matrices (also known as their Hadamard product), and EXP(.) denotes the entrywise 
exponential function acting on a matrix. 

Definition 2. Given a (dxd) complex Hadamard matrix H , an affine family of Hadamard 
matrices stemming from it is given by 



where 1Z is a subspace of all real (d x d) matrices with zeros in the first row and column. 

A family of Hadamard matrices is called non-affine if it cannot be written in this form. 

The defect d(H) of a Hadamard matrix H, defined in )19| . provides an upper bound 
on the dimensionality of any set of Hadamard matrices stemming from H. If a dephased 
Hadamard matrix has a defect of zero then the matrix is called isolated, expressing the fact 
that all complex Hadamard matrices in a neighbourhood of H are equivalent. To calculate 
the defect of H, we multiply all elements of the Hadamard matrix with independent phase 
factors, i.e. Hij — > e aij Hij for i, j = 2 . . . d, and solve the set of equations, to first order, 
which are imposed by the unitarity condition (see [27] for an explicit example) . For all but 
the smallest dimensions d or special cases, it seems imperative to use a computer program 
in order to determine the defect. 

2.2 Known constructions in composite dimensions 

There are many known constructions of complex Hadamard matrices (cf. |20| ) some of 
which apply only to specific dimensions. We briefly review two constructions of (affine) 
complex Hadamard matrices based on the tensor product of smaller Hadamard matrices. 
The two constructions described below are summarised in |28| . 

Theorem 1. (Hosoya-Suzuki [13]) Let M\,M 2 , ...,M V be k x k, N 1} N 2 , ■ ■ ■ , N k be v x v 
complex Hadamard matrices. Then the generalised tensor product matrix, denoted by Q = 
(Mi, M 2 , ■ ■ ■ , M v ) (g) (Ni,N 2 , ■ ■ ■ ,Nk), whose (i,j)th block is given by the matrix Qij = 
diag([M\]ij, [M 2 ]ij, ■ ■ ■ , [M v ]ij)Nj , is a complex Hadamard matrix of order vk. 

By using a simpler version of this tensor product structure, these matrices can be param- 
eterised, i.e. embedded in larger families of Hadamard matrices. 

Corollary 1. (Di^a [14]) Let M = (mij) be a k x k and N±, N 2 , . . . , N k be v x v dephased 
complex Hadamard matrices with m and n\,n 2 , . . . ,n k free parameters, respectively. Then 
the block matrix 



H(R) = {Ho EXP(i -R):R€TZ} 



(2) 





Q 
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with blocks Qij = rriijNj, is a complex Hadamard matrix of order vk with m + X^i=l n i + 
(k — l)(v — 1) free parameters. 

Any matrix which can be derived from Corollary [T] is called a Dif,a-type complex Hadamard 
matrix. 

2.3 Butson-type complex Hadamard matrices 

We now recall a special class of Hadamard matrices called Butson-Hadamard matrices. 

Definition 3. A complex Hadamard matrix of order d is a Butson-Hadamard matrix 
BH(d,r) if its elements are r th roots of unity, apart from a factor \j\fd. 

Note that some authors write BH(r,d) instead of BH(d,r). It is straightforward to show 
that a Hadamard matrix is (equivalent to one) of Butson-type BH(d,r): once dephased, 
all its matrix elements must be r th roots of unity. 

The simplest examples of Butson-type matrices occur when r = 2; in this case the matrices 
BH(d, 2) are the set of (dxd) real Hadamard matrices. The existence of BH(d, r) matrices 
for arbitrary values of d and r is still an open problem; it remains unknown, for example, if 
real Hadamard matrices of the form B(4n, 2) exist for all integers n. A summary of existing 
Butson-Hadamard matrices with fourth and sixth roots of unity can be found in [28j and 
a catalogue of known BH(d, r) matrices for d < 16 is given in |20] , There are also several 
existence theorems for BH(d, r) matrices, e.g. 

Theorem 2. (Butson |29| ) When p is prime, a BH(2p,p) matrix can be constructed. 

For p = 3, the matrix BH(6,3) turns out to be the isolated matrix Sq, which was also 
found independently in |30[I31|, We have derived the matrices Bio £ BH (10, 5) and 
B14 G BH(14:, 7) following Butson's method (see Appendix|Aj since they will be important 
in the present context and seem to be unavailable in the literature. 

Butson-Hadamard matrices will appear in Sections |4] and [5j where we will derive some 
previously unknown examples of complex Hadamard matrices of orders up to 91 based on 
a different construction technique. Most of these examples cannot be constructed from 
Theorem [2] 

3 Complex Hadamard matrices from pairs of MU product 
bases 

The following theorem shows how to construct a complex Hadamard matrix of order d = pq 
where p and q are both prime, from sets of smaller MU bases. 
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Theorem 3. Suppose that Kq, . . . , K p _\ and Lq, . . . , L p ^\ are unitary matrices of order 
q such that KjnL n are complex Hadamard matrices for all m, n = 0, . . . ,p — 1, i.e. K m 
is MU to L n , and let ctij/y/p be the (i,j)th element of a complex Hadamard matrix M of 
order p, with \ \ = 1 . Then the block matrix H pq given by 



H 



pq 



1 



Ol 2x K[Lq 



ct\ 2 K\hx 
a 2 2K\Li 



"31-^2^0 a 32 KlLi 



ot?, p K 2 Lp-i 



is a complex Hadamard matrix of order pq. 



(4) 



The theorem follows easily from factorising the matrix H pq such that H pq = b\b 2 , where 



/ K 
K x 



\ 



and 



Bo 



V 

/ "11-^0 

a 2 iL 





a\ 2 Li 
a 22 Li 



K 



(5) 



p-1 J 



ctipLp-i ^ 
ct 2 pL v ~\ 



(6) 



\ ctpiL a p2 Li 

The column vectors of the unitary matrices B\ and B 2 form a pair of MU bases since 
the block matrices K m are MU to L n , i.e. KmL n are complex Hadamard matrices for 
all m,n = 0,...,p— 1. Thus, by mapping B\ to the identity matrix using the unitary 
transformation B\ 1 the matrix B 2 is simultaneously mapped to B\B 2 
complex Hadamard matrix. This completes the proof of Theorem [3] 



Hpq, which is a 



In fact, the matrices B\ and B 2 correspond to a pair of MU product bases where the 
columns of B\ and B 2 form the vectors of each basis. We can write the pair of matrices 
B\ and B 2 as the orthonormal bases 



10, 



/C , \l z ) ®/Ci, 



I(P-1), 



and 



e 2 = {|0 a )®£ , |la>®A, |(p-l)o>®>Cp_i} 



(7) 



(8) 



respectively, where \m z ) ®/C m denotes the tensor product of a state \m z ) from the standard 
basis of C p with all states from a basis JC m of the space C q corresponding to the matrix 
K m . Similarly, \n a )<S>C n is defined such that \n a ) is a state in C p corresponding to the nth 
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column vector of the {p x p) Hadamard matrix M, and with C n a basis of C q corresponding 
to the matrix L n . Thus, by mapping 8\ to the standard basis, the columns of the second 
basis 82 form the complex Hadamard matrix H pq . 

To simplify the matrix H pq we perform equivalence transformations on the pair of MU 
bases {81,82} such that /Co and Co are mapped to the standard and Fourier basis of C q 
respectively, and the orthonormal basis {|0 a ), . . . , — l) a )} is mapped to the Fourier basis 
of C p , i.e. Kq = I q , Lq = F q and M = F p , with I q the {q x q) identity matrix, and F p , F q , 
being the Fourier matrices of order p and q respectively. Since B\ is MU to B2, the set 
{Iq, K\, K2, ■ ■ ■ , Kp-i} is MU to {F q , L\,..., L p _i}, and as a consequence, L\, . . . , L v ^\ 
are complex Hadamard matrices. We will continue to use the simplification M = F p , 
Kq = Iq and Lq = F q throughout. 

In the trivial case of K\ = . . . = K p ^\ = I q , one can choose each matrix L n to be a 
(q — l)-parameter family DF q where D = diag(l,e a " , . . . ,e a i- 1 ) for each n > 0. In this 
case, H P q is a (p — l)(q — l)-parameter family of complex Hadamard matrices. 

In the following section we will show that for certain choices of the (q x q) matrices 
K\, . . . , K p _i,L\, . . . , Lp-i, the matrix H pq given in Theorem [3] supplies new examples 
of Hadamard matrices. Most of the matrices we find will be isolated, which is sufficient 
to confirm that Theorem [3] produces matrices not of Di^a-type: every Di^a-type matrix 
is embedded within a family depending on at least (k — l)(v — 1) free parameters, with 
k,v > 1. 

4 Examples: d < 15 

We will now use the construction given in Theorem |3]to find complex Hadamard matrices of 
composite dimensions d = pq, with prime numbers p < q. In this section, we limit ourselves 
to matrices of order d < 15 with p < q. Larger dimensions and possible generalisations of 
the construction will be considered briefly in the following section. 

4.1 Dimension four 

In dimension four, all inequivalent complex Hadamard matrices are given by the one- 
parameter family F±(a), a £ [0,7r] j9[[T0]. We re-derive this family from Theorem [3] using 
the block matrix, 





where 




(10) 
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is the (2 x 2) Fourier matrix and L\,K\ are specific unitary matrices of order two: to apply 
Theorem [3] it is necessary that the set {I2,K\} is MU to {i*2,£i}. 

If K\ is chosen as the identity, then L\ can take the form 



where the column vectors of L\ are indeed MU to the standard basis; an overall phase factor 
has been removed using equivalence transformations. Thus, H4 turns into a one-parameter 
family of complex Hadamard matrices, 



By permuting rows it is easily shown that H^(a) is equivalent to the one-parameter Fourier 
family F±{a). One can exchange K\ with L\ but the resulting family is still equivalent to 
F±{a); no other choices are possible. Note that in the four-dimensional case, F±{a) is 
equivalent to the transposed Fourier family (F±(a)) T , a relation that does not always hold 
for larger composite dimensions. 

4.2 Dimension six 

In a recent paper (23] , pairs of MU product bases of the form B\ = {|0) (g) ATo, |1) ®K,\\ 
and Bi = {|0 X ) ® Cq, \1 x ) ® £1} were shown to give rise to the transposed Fourier family 
of complex Hadamard matrices [19] and the isolated matrix Sq E BH(6,3) [29H31] , Here 
we rederive this result on the basis of Theorem [3] where we start with (6 x 6) matrix 



three. The bases B\ and B2 will be MU if the pair {13,^1} is MU to the pair {F3, L\}. A 
proof given in |22] limits the possible choices for the matrices K\ and L\ to just three: (i) 
K\ = I3; (ii) L\ = F3; (iii) all four matrices are pairwise MU. 

If K\ = I3, the most general set of matrices satisfying the MU conditions is the two- 
parameter set, 




(11) 




(12) 




(13) 





(15) 
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Thus, Hq becomes a two-parameter family of complex Hadamard matrices which is equiv- 
alent to the transposed Fourier family (F^) T . 

If all four matrices I 3 ,F 3 ,Ki and L\ are MU then K\ and L\ are given by 



Kl = Hy = 7 = 3 



( 1 1 1 



1 

00 co 1 and L\ = H w = —=- 

12 

CO 1 00 



( 1 1 1 

(J 2 1 w I , 

V 



oo 2 00 1 



or vice versa. Here, {I 3 , F 3 , H y , H w } is the complete set of MU bases in C 3 . The complex 
Hadamard matrix Hq in (|13p associated with K\ = H y and L\ = H w is equivalent to Sq, 
the only known isolated complex Hadamard matrix of order six. 

Thus, we have indeed constructed the transposed Fourier family of complex Hadamard 
matrices and the isolated matrix Sq from Theorem [3] There is a multitude of additional 
Hadamard matrices of order six, including three- and four-parameter families |11|I12]. but 
none of these can be derived from Theorem [3] 

4.3 Dimension nine 

The online catalogue j2U] lists three types of complex Hadamard matrices of order nine: 
the four-parameter Fourier family -Fg 4 \ the isolated matrix Ng [32], and the matrix Bg 



which has a defect of two |32] . Three Butson-Hadamard matrices are among them, namely 
F 3 ® F 3 G BH(9, 3), F 9 E BH(9, 9) and B 9 G BH(9, 10). Theorem [3] allows us to identify 
an additional isolated Butson-Hadamard matrix of the form BH{9,6). 

The matrix in Eq. (|3|), for d = 9, has the structure 



Ha — — — 

V3 



' F3 L\ L2 

k\f 3 ooK\l x oo 2 k\l 2 I , 
\ K ] 2 F 3 oo 2 K\L x ooK ] 2 L 2 



where 00 = e 27 ™/ 3 is a third root of unity, and the matrices I3 and F 3 are the (3 x 3) identity 
and Fourier matrix, respectively. The (3 x 3) matrices K\, K 2 , L\ and L 2 must be chosen 
such that the set {/ 3 , I<i,K 2 } is MU to {F 3 , Li, L 2 }. 

In the six-dimensional case (cf. Sec 14 . 2 [) . the choice of pairs {I,K\} and {F 3 ,Li} was 
limited to either two matrices having identical columns (up to column permutations) or all 
four matrices being MU. Similarly, the choices for the triples {I 3 , K\, K 2 } and {F 3 , Li,L 2 } 
is restricted to the following two possibilities: (i) three matrices within one triple are 
identical (up to column permutations); (ii) two matrices are identical in each triple and all 
four MU bases {I 3 , F 3 , H y , H w } are used. 

If all matrices in one triple are the same, i.e. K\ = K 2 = I 3 , then Ha is equivalent to the 
transposed Fourier family (Fg^) T of order 9, which depends on four real parameters. 
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Now suppose that K\ = I 3 and L\ = F 3 . The only remaining choice for matrices K2 and 
L2, (if K2 7^ ^3), that satisfy the MU conditions is K2 = H y and L2 = H w (or vice versa), 
with H y and H w defined in Eq. (|16p . Denoting the resulting matrix by Sg, we find 

/ 



So 



1 

71 



F 3 F 3 H w 
F 3 loF 3 lo 2 R w 
\ HlF 3 Lo 2 H ] y F 3 u)HlH u 



(18) 



After dephasing, the matrix Sg takes the form 



So 



/ 1 


1 


1 


1 


1 


1 


1 


1 


1 


\ 


1 


LO 


LO 2 


1 


LO 


LO 2 


LO 2 


1 


LO 




1 


LO 2 


LO 


1 


LO 2 


CO 


CO 2 


LO 


1 




1 


1 


1 


LO 


LO 


LO 


LO 2 


LO 2 


LO 2 




1 


LO 


LO 2 


LO 


LO 2 


1 


LO 


LO 2 


1 




1 


lo 2 


LO 


LO 


1 


to 2 


LO 


1 


LO 2 




1 


LO 


LO 


LO 2 


1 


1 


— LO 


-1 


-1 




1 


lo 2 


1 


LO 2 


LO 


LO 2 


-LO 2 


-1 


-LO 2 




V 1 


1 


CO 2 


LO 2 


oo 2 


LO 


-LO 2 


-LO 2 


-1 


/ 



(19) 



where lo = e 2m ^ is a third root of unity. Due to the negative signs in the bottom right 
block, Sg is a Butson-Hadamard matrix containing sixth roots of unity, i.e. Sg G BH{9, 6). 

We find the defect of this matrix to be d(Sg) = implying that Sg is isolated. 
Proposition 1. The matrix Sg is inequivalent to Fg , Bg and Ng. 

Since Fg 6 ^ and Bg have non-zero defects and JVg contains only tenth roots of unity, it is 
clear that Sg is inequivalent to any known complex Hadamard matrix in dimension d = 9. 
The only other isolated complex Hadamard matrix known for d = 9, i.e. Ng, was found by 
a numerical search in |32] , As far as we are aware, the matrix Sg has not been published 
previously. 



4.4 Dimension ten 

The known Hadamard matrices in dimension ten include the Fourier family F^ and its 
transpose (F^) T , the family found by Dit^a |17] and a family stemming from 
D\q [33j. There is also an isolated matrix N\qa, a family N^ 0B found originally in |32] 
and parameterised in |34] . and G^j |34| . Furthermore, there are the Butson-Hadamard 
matrices Xiq E BH(W, 5) [35] and W' G BH(10, 6) [36J. Within the continuous families, 
several Butson-type matrices exist: D\q G BH(10, 4), F2 <8> F5 ~ F10 € i?i/(10,10), and 
those contained in , e.g. h[q\ d^ € BH (10, 6) given in |17| . 
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In the following, we construct a complex Hadamard matrix of Butson-type based on the 
block matrix 



Hi n — —= 



F 5 
•t 



V2 \ K[F 6 —K\Li 



(20) 



it is necessary that the pair of (5 x 5) matrices {I5, K\} is MU to the pair of (5 x 5) matrices 
{F^,L\}. Here, 1$ is the identity matrix and F§ the Fourier matrix, 



F R = 



V5 



( 1 


1 


1 


1 


1 


\ 


1 


UJ 


c^ 2 


a; 3 


uj [ 




1 


uj 2 


UJ 1 


UJ 


UJ 3 




1 


UJ 3 


UJ 


uj 4 






V 1 


uj 4 


UJ 3 


UJ 2 


UJ 


/ 



(21) 



with fifth root of unity. Assuming that K\ and L\ are not identical to I5 and 

F$, respectively, one choice is to require that the matrices within the set {1$, F$, K\, L\] 
are pair wise MU. 

In dimension five, the complete set of six MU bases can be written as 

{-^5,-^5) Hi,H 2 , H 3 , H4}, 
where Hi are the complex Hadamard matrices of order five given by 

H 1 = DF 5 , H 2 = D 2 F 5 , H 3 = D 3 F 5 , H A 
and with a diagonal matrix, 

D = diag(l, uj, uj 4 , uj 4 , uj). 



D 4 F 5 , 



(22) 



(23) 



(24) 



This characterisation of the complete set of MU bases is based on a construction in [TH 
By choosing K\ = H3 and L\ = Ha, the matrix H\q becomes 



Sw — —= 



Ft 



Ha 



y/2 \ HlF 5 -HIH4 



(25) 



which has the dephased form 



S 
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( 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 






1 


UJ 


UJ 2 


UJ 3 


UJ 4 


UJ 4 


1 


u 


UJ 2 


UJ 3 






1 


UJ 2 


UJ 4 


UJ 


UJ 3 


■UJ 


UJ 3 


1 


UJ 2 


UJ 4 






1 


UJ 3 


UJ 


UJ 1 


UJ 2 


UJ 


UJ 4 


c^ 2 


1 


UJ 3 






1 


UJ 4 


UJ 3 


uj 2 


UJ 


UJ 4 


uj 3 


uj 2 


UJ 


1 






1 


UJ 3 


uj 2 


uj 2 


uj 3 


1 


UJ 


UJ 4 


uj 1 


UJ 






1 


UJ 2 


1 


UJ 1 


UJ 4 


UJ 3 


UJ 2 


UJ 3 


UJ 


UJ 






1 


UJ 


uj 3 


UJ 


1 


UJ 2 


UJ 4 


UJ 3 


UJ 4 


UJ 2 






1 


1 


UJ 


UJ 3 


UJ 


UJ 2 


UJ 2 


UJ 4 


UJ 3 


UJ 4 




V 


1 


UJ 1 


UJ 1 


1 


UJ 2 


UJ 3 


UJ 


UJ 


UJ 3 


UJ 2 


) 



(26) 
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where uj = e 27 ™/ 5 is a fifth root of unity. 

Proposition 2. The matrix S\q is inequivalent to F±q , (F^) T , , I^q , A^oA; -^xob> 
G$ and W'. 

The matrix Siq is found to be isolated and contains only fifth roots of unity, therefore it is 
inequivalent to any of the complex Hadamard matrices listed in the proposition. However, 
we have not been able to show whether it is equivalent (or not) to the isolated Butson-type 
matrix Xiq £ BH(10, 5) or the matrix B\q given in Appendix lAl 

A different choice of the matrices K\ and L\, for example, K\ = Hi and L\ = H 2 , leads 
to the block matrix 



5' -J- 



F 5 
rt 



H 2 

t; 



H{F 5 -H[H 2 



(27) 



which in dephased form is written as 



10 





1 


1 


1 


1 


1 


1 


1 


1 


1 


1 






1 


Co' 


Co- 2 


Co' 3 




UJ 2 


UJ 3 


OJ 4 


1 


Co' 






1 


UJ 2 


UJ 4 


Co' 


UJ 3 


UJ 3 


1 


UJ 2 


Co' 1 


Co' 






1 


UJ 3 


u) 


OJ 4 


UJ 2 


to 3 


UJ 


Co' 1 


Co' 2 


1 






1 


Co' 1 


uj 3 


uj 2 


UJ 


UJ 2 


UJ 


1 


Co- 4 


uj 3 






1 


Co' 1 


UJ 


UJ 


uj 4 


-1 


— UJ 


-UJ 4 


-uj 4 


— UJ 






1 


UJ 


1 


Co' 2 


UJ 2 


-UJ 2 


— UJ 


-UJ 2 


-1 


-1 






1 


UJ 3 


Co' 1 


uj 3 


1 


-UJ 3 


-1 


-UJ 4 


-1 


-Co- 3 






1 


1 


UJ 3 


Co' 1 


UJ 3 


-UJ 3 


-u; 3 


-1 


-uj 4 


-1 




{ 


1 


Co- 2 


UJ 2 


1 


UJ 


-UJ 2 


-1 


-1 


-UJ 2 


— Co! 


J 



(28) 



with fifth root of unity. This matrix is a member of the family BH(W, 10) 

and, with a defect d(S' w ) = 8, the maximum dimension of any smooth manifold stemming 
from S' w will be eight. Several other matrices of the form BH(W, 10) exist but we are not 
able to determine whether S[ is equivalent to any of them. 

Any choice of K\ and L\ from the set of MU bases {-F5, Hi, H 2 , H3, H4} will result in 
a Butson-type matrix of the form BH(10,5) or BH(10, 10). It would be interesting to 
see if the various combinations of Ki and Li result in further new inequivalent complex 
Hadamard matrices. 



4.5 Dimension fourteen 

The known complex Hadamard matrices of order fourteen are the six-parameter Fourier 
family F$ and its transpose (F^) T , the family Dy± found in [33], and a set of isolated 



12 



matrices L^ x for X 



A,B,C, . . . ,N found in [34] . In addition, there are several Di^a- 
type matrices, listed in [20], obtained from Diva's method of Corollary [1] using the Fourier 
matrix F2 and any Hadamard matrix of order seven. 



The matrix we construct from Theorem [3] consists of four blocks, 

1 



H 



14 



F 7 
•t, 



-t 



V2 \ K\F 7 -K[L x 
where the Fourier matrix of order seven is given by 



(29) 





1 


1 


1 


1 


1 


1 


1 


\ 




1 


OJ 


OJ 2 




oj 4 


OJ 5 








1 


OJ 2 


oj 4 


OJ 6 


OJ 


OJ 3 








1 


OJ 3 


OJ 6 


OJ 2 


OJ 5 


OJ 


w 4 






1 


OJ 4 


oj 


oj 5 


LO 2 


OJ 6 


a; 3 






1 


oj 5 




OJ 


OJ 6 


Co' 1 


a: 2 




I 


1 


OJ 6 


OJ 5 


oj 4 


OJ 3 


UJ 2 


OJ 


J 



(30) 



with oj = e 2nl / 7 , and the (7 x 7) matrices K\ and L\ are chosen from the complete set of 
eight MU bases of the space C 7 . 



Let us denote the complete set of MU bases by F7, H\, H2, H^,!!^, H^,Hq} where 



and 



Hi = D 3 F 7 



D = diag(l, 1, oj, oj 3 , oj g , oj 3 , oj). 



(31) 
(32) 



with j = 1, . . . , 6. The diagonal D is based on the construction of a complete sets of MU 
bases in prime dimensions presented in |37j . By choosing K\ = H\ and L\ = H2, we find 
a complex Hadamard matrix which, after dephasing, reads explicitly 



Si 



( 1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 




oo 2 


OJ 3 


OJ 4 


LJ 5 


LO 6 


1 


OJ 




OJ 3 


OJ 4 


OJ 5 


u e 


1 


^ 


UJ 1 


OJ 6 


LO 


OJ 3 


0/' 


2 

OJ 


OJ 4 


OJ 6 


LO 


OJ 3 


OJ 5 


1 


1 


OJ 3 


0J ( ' 


OJ 2 


OJ 5 


LO 


00+ 


OJ 6 


oo 2 


OJ 5 


LO 


OJ 4 


1 


OJ 3 


1 


UJ 4 


to 


OJ 5 


2 

OJ 


OJ 6 


OJ 3 


OJ 5 


2 

OJ 


OJ 6 


OJ 3 


1 


OJ 4 


LO 


1 


UJ 5 


UJ 3 


00 


00" 


OJ 4 


oo 2 


OJ 6 


OJ 4 


OJ 2 


1 


OJ 5 


OJ 3 


OJ 


1 


OJ 6 


00" 


LO 4 


OJ 3 


^ 


LO 


OJ 2 


OJ 


1 


OJ 6 


OJ 5 


OJ 4 


OJ 3 


1 


UJ 


UJ 3 


OJ 6 


OJ 3 


LO 


1 


OJ 5 


OJ 5 


OJ 4 


LO 2 


OJ 6 


^ 2 


OJ 4 


1 


1 


LO 


OJ 3 


LO<> 


OJ 3 


OJ 


OJ 4 


OJ 5 


OJ 5 


OJ 4 


OJ 2 


OJ 6 


^ 2 


1 


OJ 6 


OJ 6 


1 


^ 2 


OJ 5 


LO 2 


OJ 


OJ 3 


OJ 4 


LO 1 


OJ 3 


LO 


OJ 5 


1 


UJ 5 


00+ 


LO 4 


OJ 5 


1 


OJ 3 


OJ 3 


OJ 6 


LO 


LO 2 


UJ 2 


LO 


OJ 6 


1 


OJ 4 


oo 2 


00 


LO 




OJ 4 


OJ 3 


1 


OJ 3 


00" 


OJ 8 


LO (l 


OJ 5 


1 


UJ 3 


1 


OJ 5 


OJ 4 


to 4 


OJ 5 


LO 


OJ 6 


OJ 3 


LO" 


LO 


2 

OJ 


LO 2 


^ 1 


00 2 


00" 


UJ 2 


1 


J' 


0J« 


OJ 4 


OJ 3 


LO 


0/' 


OJ 


OJ 3 


OJ 4 



(33) 
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with u = e 2 ™/ 7 . This is an isolated Butson-type complex Hadamard matrix of the form 
BH(U,7). 

Proposition 3. The matrix S14 is inequivalent to F^ , (F^) T , Dy± , L^ x and every 



Difa-type matrix given in \2ty . 

Since Si 4 and L^ x contain different roots of unity they are inequivalent. All other known 
matrices listed in the proposition are contained in families of complex Hadamard matrices, 
thus, they are inequivalent to S14. However, it is not known whether S14 is equivalent to 
the Butson-Hadamard matrix B14 given in Appendix lAl 

We can construct further complex Hadamard matrices by choosing different MU bases for 
K\ and L\, e.g. if K\ = H\ and L\ = H4, the resulting matrix in dephased form is 

\ 



S' u — 



1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


UJ 


^ 2 


OJ 3 


UJ 1 


UJ 5 


oo 6 


1 


to 


UJ 2 


UJ 3 


UJ 4 


to 5 


UJ 6 


1 


UJ 2 


UJ 4 


UJ 6 




to 3 


OJ 5 


OJ 4 


OJ 6 


'jj 


u 3 


to 5 


1 


UJ 2 


1 


UJ 3 


UJ 6 


uj 2 


to 5 


OJ 


UJ 4 


UJ 5 


to 


UJ 4 


1 


to 3 


u 6 


uJ 2 


1 


to 4 


to 


UJ 5 


UJ 2 


UJ 6 


CO 3 


OJ 3 


1 


to 4 


UJ 


UJ 5 


UJ 2 


to 6 


1 


LJ 5 


UJ 3 


UJ 


OJ 6 


LJ 4 


CO 2 


UJ 5 


a; 3 


to 


co e 


UJ 4 


UJ 2 


1 


1 


u e 




UJ 4 


to 3 


LJ 2 


to 


CO 4 


a; 3 


UJ 2 


UJ 


1 


UJ 6 


UJ 5 


1 


to 


CO 3 


UJ 6 


to 3 


UJ 


1 


-UJ 3 


— UJ 


—UJ 


-UJ 3 


-1 


-UJ 6 


-I 


1 


1 


to 


UJ 3 


LO 6 


UJ 3 


to 


-1 


-UJ 3 


—UJ 


—UJ 


-UJ 3 


-1 


-UJ 6 


1 


CO 6 


UJ 6 


1 


UJ 2 


CO 5 


UJ 2 


-UJ 5 


-iu e 


-UJ 2 


-I 


-1 


-UJ 2 


-UJ 6 


1 


LJ 5 


CO 4 


UJ 4 


CO 5 


1 


CO 3 


-UJ 4 


-UJ 3 


-UJ 4 


-I 


-UJ 5 


-UJ 5 


-1 


1 


LJ 4 


uj 2 


UJ 


UJ 


LJ 2 


UJ 4 


-UJ 4 


— UJ 


-I 


—UJ 


-to 4 


-UJ 2 


-UJ 2 


1 


UJ 3 


1 


OJ 5 


to 4 


OJ 4 


OJ 5 


-UJ 5 


-I 


-UJ 4 


-UJ 3 


-UJ 4 


-I 


-UJ 5 


1 


u 2 


UJ 5 


uj 2 


1 


LJ 6 


CO 6 


-1 


-I 


-UJ 2 


-UJ 6 


-UJ 5 


-UJ 6 


-UJ 2 



J 

(34) 

This is a Butson-Hadamard matrix of the form BH^IA, 14), and has a 
defect of 12. It is unknown if S 14 is equivalent to a BH(14, 14) matrix contained within 
an existing family of Hadamard matrices. 



with co 



4.6 Dimension fifteen 



The only known complex Hadamard matrices of order fifteen are the eight-parameter 



<E £#(15,15), and 

the transposed Fourier family (F^) T . To construct a new (15 x 15) complex Hadamard 



Fourier family F 15 , stemming from the Fourier matrix F15 
F01 

matrix by means of Theorem [3j we use the block matrix 



H 



15 



1 

71 



/ F 5 Li L 2 
K\F s uK\Li a 2 K[L 2 



(35) 
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where a = e 27 ™/ 3 is a third root of unity, ^5, F5, K\, Ki, L\ and hi are (5 x 5) matrices 
including the identity matrix I5 and the Fourier matrix F5 defined in Eq. (121 p . The set 
{I 5 ,K 1 ,K 2 } is MU to the set {F 5 ,hi,h 2 }. 

If we use the complete set of six MU bases the space C 5 , {I5, F5, Hi, Hi, H3, H4}, corre- 
sponding to Ki = Hi, K<i = Hi, hi = H3 and hi = H±, as defined in Eq. (|23p . the 
resulting complex Hadamard matrix becomes 



' ^5 ^3 Hi 

h\f 5 aH\H 3 a 2 H\H A 

Apart from a factor l/\/T5, its dephased form reads explicitly 



1 

7s 



(36) 



V 


15, its 


1 


1 


u> 12 








9 


^ 6 




u> B 








u> 9 






u> 12 


-u, 5 






u> B 


-oj 2 






1 




^ 


-co 8 








-u 2 








_ w 10 


u, 3 








1 


_ w 13 




_ w 13 


, ,9 





(37) 

where u; = e 27ri ' 15 is now a fifteenth root of unity. Since all matrix elements can be 
written in terms of 30th roots of unity, S15 is an example of a Butson-Hadamard matrix 
BH (15, 30). The vanishing defect of this matrix, i.e. d(Sis) = 0, implies that S15 is 
isolated. This property excludes S15 from being a member of either of the affine families 
F® and (Fiff. 

Proposition 4. The matrix S15 is inequivalent to F^' and (F^) T . 



One could produce additional complex Hadamard matrices by choosing different combina- 
tions of MU bases from the complete set of six, such as Ki = 1$ or hi = F3. It is likely 
that various inequivalent matrices will result from these choices. 



5 Examples: dimensions d > 15, and further generalisations 

The construction of the matrices Sq, Sg and S15 has a common feature: in each case, the 
matrices Kq, . . . , K p _i,ho, . . . , L„_i— used to construct the blocks of the Hadamard matrix 
in Theorem |3]-include a complete set of (q + 1) MU bases of the space C q . For the cases 
d = 6 and d = 15, i.e. q = 2p — 1, complete sets of MU bases in dimension three and 
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five are used, respectively, resulting in the isolated matrices Sq and S15. Furthermore, in 
the case d = 9, where q < 2p — 1, a complete set of four MU bases in dimension three is 
used, and again we find an isolated Hadamard matrix Sg. One might expect additional 
isolated complex Hadamard matrices to emerge for larger composite dimensions whenever 
its factors are related by q < 2p — 1. 

We have been able to confirm this property for all primes p, q, with pq < 100 and q < 2p—l, 
excluding the case d = 4. The first three examples are covered by 5*g, Sg and S15 which 
we already know are isolated. The remaining five matrices S25, S35, S49, S77 and Sgi, also 
turn out to be isolated. We construct these matrices as follows: 

• S25 is derived from the complete set of six MU bases I5, F$ and Hj = D 3 F^, j = 1 ... 4, 
where D = diag(l, u), uj 4 , uj 4 , uj) and u> = e 27 ™/ 5 . The matrices K n , for n = . . . 4, are 
chosen as 1$, I5, I5, Hi, H 2 , respectively, and L n as F§, F§, F§, H3, H^, respectively. 

• S35 uses the complete set of eight MU bases I 7 , F 7 and Hj = D J Fj, j = 1 ... 6, where 
D = diag(l, 1, uj, uj 3 , uj 6 , uj 3 , uj) and uj = e 2nl ^ 7 . The matrices K n , for n = ... 4, are chosen 
as Hi, H2, H3, and L n as -F7, -F7, H^, H^, Hq, respectively. 

• S49 is constructed from the same complete set of MU bases used for S35, and we choose 
K n as I 7 ,I 7 ,I 7 ,I 7 ,Hi,H 2 ,H 3 , and L n as F 7 , F 7 , F 7 , F 7 , H 3 , H A , H 5 , with n = ... 6, re- 
spectively. 

• S77 uses a complete set of 12 MU bases given by In, F\\ and Hj, j = 1 ... 10 with 
D = diag(l, 1,oj, uj 3 ,uj e , uj 10 , uj 4 , uj 10 , uj 6 ,uj 3 ,uj) and uj = e 27 ™/ 11 . The matrices K n , for n = 
0...6, are chosen as In, In, Hi, H 2 , H 3 , i? 4 , H 5 and L n as Fn, F n , H 6 , H 7 , H$, H 9 , H w , 
respectively. 

• 5gi is based on the complete set of 14 MU bases in C 13 , i.e. the identity /13, the Fourier 
matrix F13, and the matrices Hj = D J Fi 3 for j = 1 ... 12, where the diagonal matrix 
is given by D = diag(l, 1, uj, uj 3 , uj 6 , uj 10 , uj 2 , uj 8 , uj 2 , uj 10 , uj 6 , uj 3 , uj) and uj = e 2m / 13 . The 
matrices K%, . . . , Kq,Li, . . . ,Lq, correspond to Hi, . . . , Hyi, respectively. 

All of these isolated matrices are of Butson-Hadamard type with S25 £ BH(25, 10), S35 £ 
B#(35,70), 5 49 G £#(49,14), S77 G BH(77, 154) and S91 G BH (91, 182). They may 
have smaller roots of unity if their matrix elements contain no entries equal to (—1). The 
matrix Sgi is similar to Sq and S15 in the sense that the prime factors of d = 91 satisfy 
the equality q = 2p — 1, meaning that each MU basis from the complete set is used exactly 
once. For the other isolated matrices, the factors satisfy the inequality q < 2p — 1, which 
implies that some MU bases are used more that once in the set Kg, . . . , L p _i. In this case, 
there are additional choices for the bases used; different combinations may lead to further 
inequivalent isolated complex Hadamard matrices. 

So far, we have applied Theorem [3] mainly to product dimensions 4 < d < 15 or when 
d = pq < 100 and q < 2p— 1. To explore whether the latter constraint on the factors p and 
q is necessary, we have constructed Hadamard matrices in all other composite dimension 
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d < 100 for p,q < 13. In each of these cases, i.e. d = 21, 22, 26, 33, 39, 55 and 65, 
we were able to identify isolated Hadamard matrices. In addition, it is also possible to 
construct Hadamard matrices with non-zero defects, simply by selecting different sets of 
MU bases for Kq, . . . , L„_i. Thus, the theorem is potentially the source of infinitely many 
new Hadamard matrices in arbitrary product dimensions. 

Interestingly, the method is not limited to dimensions of the form d = pq: if the numbers p 
and q are composite, it is likely that additional, possibly inequivalent complex Hadamard 
matrices can be constructed which relate to different factorisations of the dimension, such 
as2x6 = 3x4 when d = 12. Furthermore, it has been shown |38| that inequivalent 
complete sets of MU bases exist for large prime powers d = p n , possibly leading to yet 
more inequivalent Hadamard matrices. One could also try to create continuous families 
of complex Hadamard matrices if sets of four or more MU bases exist which contain free 
parameters after dephasing. 

Finally, another generalisation of Theorem [3] can be achieved as follows. The Hadamard 
matrices we have constructed are derived from product bases which tensor each vector in an 
orthonormal basis of C p with an orthonormal basis of C q (cf. Eqs. (|7|8p ). However, other 
types of product bases exist; for example, one could take vectors from different orthonormal 
bases in C p and tensor them with vectors from one basis in C q . The classification of all 
product bases in the space C 2 (8> C 3 , up to local equivalence transformations, contains a 
number of examples of these so-called indirect product bases [22j. Thus, alternative block 
structures may be allowed in Theorem [3l potentially leading to other Hadamard matrices. 

6 Summary and outlook 

The main results of this paper are (i) a new general construction of complex Hadamard 
matrices in composite dimensions d = pq (p, q prime) described in Theorem [3l and (ii) 
the explicit derivation of various new complex Hadamard matrices as a consequence of this 
theorem. The construction, detailed in Theorem |3l relies on the simple idea that a suitable 
unitary transformation maps a pair of MU product bases to its standard form in which 
the vectors of one basis become the columns of a complex Hadamard matrix. It becomes 
possible to systematically construct new Hadamard matrices, while previous examples have 
been found by trial and error [30] or from numerical methods |32| . 

To illustrate the approach we first derive some known results in low dimensions. In par- 
ticular, we find the complete family of complex Hadamard matrices when d = 4, and 
in dimension six we find the isolated matrix Sq. We then proceed to higher dimensions, 
obtaining isolated Hadamard matrices of order 9, 10, 14 and 15. Two of these are new 
isolated Butson-type Hadamard matrices, namely Sg G BH{9,6) and S15 G -B-ff (15, 30), 
while S10 G BH(10, 5) and S14 G -6(14, 7) are shown to be inequivalent to nearly all known 
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Hadamard matrices of their order. However, we cannot exclude the equivalence of Siq to 
i?io or -XiOi an< i 5'i4 to S14. 

In dimensions d = 10 and d = 14, there is some flexibility in selecting suitable subsets of MU 
bases when applying Theorem [3] This enables us to construct two non-isolated matrices 
S' 1Q and S' 14 , with defects equal to 8 and 12, respectively. Further research is needed to 
understand which choices of MU bases will lead to inequivalent Hadamard matrices. 

Whenever the factors in the product dimension d = pq are related by q < 2p — 1, the set 
Kq, . . . , Kp-i, Lq, . . ., Lp-i given in Theorem [3] can accommodate a complete set of MU 
bases for the space C q . We speculate that in these cases, with the exception of dimension 
four, Theorem [3] will always give rise to an isolated Hadamard matrix. This expectation 
is confirmed for all matrices of order d = pq < 100 that satisfy q < 2p — 1. In these 
cases, the matrices Sq, Sg, and Si 5, as well as S25, S35, S49, S77 and S91 all turn out to be 
isolated, and they include the largest known examples of isolated Hadamard matrices (as 
far as we know). What is more, we are also able to generate isolated and (non-isolated) 
Hadamard matrices for dimensions d = 21, 22, 26, 33, 39, 55 and 65, giving rise to a total 
of 16 isolated complex Hadamard matrices. Twelve of them are new, while the remaining 
four, namely S\o, S14, S22, £26 £ BH(2p,p), may be equivalent to matrices resulting from 
Butson's construction. 

Throughout this paper we have limited our search to only Butson-type Hadamard matrices. 
However, the general method given in Theorem [3] covers a much wider class of complex 
Hadamard matrices. We expect that many other examples of more general Hadamard 
matrices can be found by extending the choice for the unitary matrices Kn and Ln- 
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A Appendix 



In this Appendix we list the two Butson-type Hadamard matrices BH(2p,p) of order 10 
and 14, which we derive from a construction method given in [29] . 

• For p = 5, the dephased matrix BH(W, 5) is given as 



B 



10 



( 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 




1 


LO 


CO 2 


oo 3 


00 1 


oo 2 


00 1 


LO 


a; 3 


1 




1 


OJ 2 


lo 4 


00 


oo 3 


oo 3 


oo 2 


00 


1 


a; 4 




1 


LO 3 


00 


OJ 4 


oo 2 


oo 3 


LO 4 


1 


a-' 


co> 2 




1 


LO 4 


LO 3 


oo 2 


00 


oo 2 


1 


LO 3 


a; 


w 4 




1 


oo 3 


LO 2 


to 2 


oo 3 


1 


00 


00 1 


cu 4 


Co' 




1 


oo 2 


1 


oo 4 


oo 4 


oo 


00 






LO 3 




1 


00 


00* 


00 


1 


oo 4 


Of' 




LO 2 


oo 2 




1 


1 


00 


oo 3 


LO 


00 1 


OJ 2 


C' 2 


00 1 


oo 3 


\ 


1 


LO 4 


00 1 


1 


oo 2 


00 


oo 3 


C' 2 


a; 3 


00 



3 27ri/5 



is a fifth root of unity. The defect of -Bio is zero. 



where oo 

• For p = 7, the dephased matrix BH(\A, 7) is 



(38) 



/ 1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 \ 


1 


LO 


9 

LO 


LO 3 


LO 4 


LO 5 


LO 6 


LO 6 


LO 2 


LO 5 


LO 


LO 4 


1 


OJ 3 


1 


oo 2 


LO 4 


L0 G 


OJ 


LO 3 


LO 5 


LO 3 


OJ 2 


00 


1 


oo" 


LO 5 


oo 4 


1 


oo :i 


LO 6 


LO 2 


LO 5 


LO 


LO 4 


LO 5 


1 


OJ 2 


LO 4 


LO 6 


lO 


LO 3 


1 


LO 4 


LO 


LO 5 


0,' 2 


LO 6 


LO 3 


LO 5 


LO 3 


LO 


LO 6 


LO 4 


LO 2 


1 


1 


LO 5 


LO 3 


00 


oo" 


LO 4 


oo 2 


LO 3 


00+ 


LO 5 


oo" 


1 


LO 


oo 2 


1 


oo" 


LO 5 


LO 4 


LO 3 


2 

LO 


LO 


LO 6 


LO 3 


1 


LO 4 


OJ 


LO 5 


2 

LO 


1 


00+ 


oo 2 


00 


LO 


LO 2 


LO 4 


1 


LO 5 


LO 6 


LO 3 


LO 3 


oo" 


OJ 5 


1 


00 


LO 3 


LO 6 


LO 3 


LO 


1 


LO 4 


oo" 


LO 4 


LO 5 


LJ 2 


LO 2 


OJ 5 


1 


oo" 


LO 4 


LO 4 


LO 5 


1 


LO 3 


^ 


LO 


LO 3 


00 


LO 2 


oo" 


oo" 


1 


oo 2 


LO 5 


LO 2 


1 


LO 8 


oo" 


00 


00+ 


LO 3 


LO 5 


LO 3 


00+ 


LO 


1 


oo" 


LO 6 


1 


LO 2 


LO 5 


^ 2 


LO 


OJ 


LO 4 


LO 3 


LO 5 


LO 3 


LO 4 


1 


LO 3 


1 


J' 


00+ 


LO 4 


LO 5 


LO 2 


oo" 


LO 6 


LO 2 


LO 


LO 3 


LO 


V 1 


1 


OJ 


LO 3 


oo" 


LO 3 


OJ 


LO 4 


oo" 


OJ 2 


LO 2 


of' 


00 + 


oo" ) 



(39) 



with oo = e 2m l 7 being a seventh root of unity. This matrix also has zero defect. 
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